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In this work we investigate the supersymmetric version of the valence bond solid (SVBS) state. In 
one dimension, the SVBS states continuously interpolate between the valence bond states for integer 
and half-integer spin chains, and they generally describe superconducting valence bond liquid states. 
Spin and superconducting correlation functions can be computed exactly for these states, and their 
correlation lengths are equal at the supersymmetric point. In higher dimensions, the wave function 
for the SVBS states can describe resonating valence bond states. The SVBS states for the spin 
models are shown to be precisely analogous to the bosonic Pfaffian states of the quantum Hall 
effect. We also give microscopic Hamiltonians for which the SVBS state is the exact ground state. 
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I. INTRODUCTION 

Quantum antiferromagnetism offers basic paradigms 
for different phases of strongly interacting quantum sys- 
tems [H, @, H|. In addition to a rich array of classically 
ordered states, including multiple sublattice Neel order 
and non-collinear states, there are several different types 
of quantum disordered states: valence bond (VB) solids, 
valence bond liquids, dimer solids, etc. By tuning various 
couplings, one can pass through quantum phase transi- 
tions which separate these states. A class of supercon- 
ductors, including the high-T c cuprates, is obtained by 
doping the Mott insulating states with quantum antifer- 
romagnetic order. In one theoretical approach, super- 
conductivity arises from doping the valence bond liquid 
state [4]. In another theoretical approach, the supercon- 
ducting state is obtained from a symmetry rotation of 
the quantum antiferromagnetic state [B| . In this work we 
construct supersymmetric extension of the valence bond 
solid state. In particular, we show that the supercon- 
ducting valence bond liquid state can be naturally ob- 
tained from the supersymmetric rotation of the valence 
bond solid state. Our results give a mathematical precise 
validation of the above-mentioned ideas. 

We investigate extensions of the valence bond solid 
states defined by Affleck, Kennedy, Lieb, and Tasaki 
(AKLT) On any lattice £, one can define a one- 

parameter family of such states, indexed by an integer 
M. The general AKLT state is written [|| 
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as a product over links (ij) of C, where S i = ^b\ 
is the local quantum spin operator, written in terms of 
Schwinger bosons, satisfying [6^,6] J = 8^6^. The 
state $(£, M)) describes an antiferromagnet where each 
site contains a single spin 5 — \zM object, where z is 
the lattice coordination number. What is special about 
these states is that the total spin J on any link can only 
take values between and J* = (z — 1)M , and all total 
spin components J = J* + 1, . . . , 25* along any link are 
absent in the AKLT wave function, because the operator 
= £fiub\ b^ v transforms as an SU(2) singlet. Thus, 
the AKLT states are annihilated by certain projection 
operators, 



spin per site is somewhere between 5 — \ and 5. 



As 



Pj(ij)\*(£,M)) =0 



(2) 



for J e [J* + 1, 25] , where J* = (z — l)M. This allows 
one to construct local Hamiltonians of the form 



2.S' 



(3) 



with Vj > 



0, which renders |<E>(£,M)) an exact, zero- 
energy ground state. In this respect, the AKLT states 
are analogous to the Laughlin wave functions in the frac- 
tional quantum Hall effect (QHE) [9j, which are also 
rendered exact eigenstates of a corresponding "truncated 
pseudopotential" Hamiltonian [j.llo|. The SU(2) AKLT 
model has been generalized by introducing g-deformcd 
SU(2) group fil l 12 L 13j|. a nd higher symmetric groups, 
such as SU(n)[!lTll El E3, SP(ra) 0, and SO(n) 

mill. 

The states we shall discuss are supersymmetric gener- 
alizations of the AKLT states, and are written as 



M, r)) = J] (e M „ blb] v + rfjf}) M \0). 



(4) 



Here, creates a fermionic hole on site i, which dis- 
places one of the bosons. The local Hilbert space thus 
accommodates two types of states: states with spin 
5 = \zM and states with spin 5 — h, and the operator 



r — > oo, each site must contain a hole, and the state is 
insulating once again. For a one-dimensional (ID) chain, 
with M = 1, there are only two possibilities: 



(5a) 
(5b) 



corresponding to spin-Peierls order. These are the two 
degenerate ground states of the well-known Majumdar- 
Ghosh Hamiltonian [U [2^, HH . In the thermodynamic 
limit, or on a ring with an even number of sites, the 
r — > oo SVBS state is the sum |$ A ) + \& B ), which has 
crystal momentum k = 0. 

On the two-dimensional square lattice, once again the 
r = state is the 5 — 2 AKLT valence bond solid. For 
r — > oo, though, rather than there being only two con- 
figurations which contribute to the SVBS wave function, 
the state is a linear combination of the resonating valence 
bond (RVB) kind, but for 5 — I. The situation is de- 
picted in Fig. [TJ The configurations which contribute to 
the SVBS state in this limit are given by dimer coverings 
of the square lattice, where each dimer corresponds to 
a hole-pair-creation operator fjfj- The quantum dimer 
model for S = h was constructed by Rokhsar and Kivel- 
son [24]]. The partition function for the classical dimer 
gas, with different fugacities for x-directed and y-directed 
dimers, was worked out by Fisher in 1961 and shown to 
take the form of a Pfaffian (25|. This connection to the 
Pfaffian is present in our work as well, and underlies re- 
cent work by one of us 0, [U [H[ on supersymmetric 
extensions of the quantum Hall problem, in which the 
Pfaffian QHE state at v — | appears as a natural limit. 
The following diagram sketches these basic connections: 

AKLT {S = 1 chain) — > Majumdar- Ghosh (S = |) 

I I 
AKLT (5 = 2 square) — > RVB (5 = §) 



Eui/ b\,,b\,. +rfjf], which creates a linear combi- T , 1 , x „, „ , \ c ■ \ 

F v i", J,J J' . . Laughlm [v = ± bosons) — > Pfaffian (v = ± fcrmions) 



ft ft 



nation of spin singlets and hole pairs on the link (ij), 
transforms as a singlet under the superalgebra OSp(l|2). 
We call these states supersymmetric valence bond solid 
(SVBS) states. Physically, the spin 5 states can be re- 
alized by 25 electrons coupled through Hund's rule cou- 
pling, and the spin 5— ^ states are obtained by removing 
one electron from the site. Thus the SVBS states de- 
scribe a doped spin chain with large on-site Hund's rule 
coupling. 

The parameter r interpolates between two limits. At 
r = 0, there are no holes, and we recover the AKLT 
state, which is an antiferromagnetic insulator. For fi- 
nite r, there is a nonzero density of nearest-neighbor hole 
pairs and the system is a superconductor. The average 



While the literature on hole motion in quantum anti- 
ferromagnets is voluminous rather little has been done to 
date to explore extensions of models of the AKLT type, 
i.e., to find wave functions at finite hole concentration 
which are exact eigenstates of local projectors (so-called 
"Klein models" [29|]). Single hole motion in the 5 = 1 
AKLT chain was discussed in Ref. [3fJ, but a differ- 
ent constraint was used: 6^ &„ + 2pf — 2. Experiments 
on hole doped AKLT spin chains have been reported in 
Ref. [3l|, where a hole divides the 5=1 AKLT chain to 
two semi- indefinite segments with 5 = \ spins at each 
edge. Interestingly, such property is shared in the SVBS 
model developed in this paper. 
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becomes 



AKLT (5 = 2) 



RVB (S = 



FIG. 1: Examples of the square lattice supersymmetric va- 
lence bond solid state with M = 1 on the square lattice. Left 
panel: r = 0, corresponding to the 5 = 2 AKLT state. Right 
panel: example configuration from the r — oo state which is 
a S = § nearest-neighbor RVB state. 



The i-J models with SU (2| 1) symmetry are known as 
the supersymmetric t-J models. The models are exactly 
solvable in one dimension [32|, [H, HH, and their corre- 
lation functions are also derived in Ref. 35]. With 1/r 2 
long-range interaction, the supersymmetric t-J models 
are still exactly solvable [36] . The models which we deal 
with possess OSp(l|2) symmetry, and their exact ground 
states are constructed even in higher dimensions as the 
case of the original AKLT models. 

The remainder of this paper is structured as follows. 
In Sec [11] we will briefly discuss the local Hilbert space 
and some preliminary aspects of the OSp(l|2) operator 
algebra, a fuller treatment of which we consign to Ap- 
pendixE] In Sec lIIII we will focus on SVBS states in one- 
dimensional systems, i.e., supersymmetric spin chains. 
Using the spin-hole coherent states developed by Auer- 
bach [371 ], we will compute various correlation functions 
in the SVBS chains. Section QV] discusses some connec- 
tions with the quantum Hall effect. In Sec|V] we derive 
a Hamiltonian with local interactions which renders our 
M = 1 SVBS chain as an exact nondegenerate ground 
state. 



II. LOCAL HILBERT SPACE AND OSp(l|2) 

In the Schwinger representation of SU(2), a spin is rep- 
resented by two bosons, with the quantum spin operator 
given by S = j&J, a b v . The total boson occupancy is 
constrained, 



a* a + b^b = p 



(6) 



where p is an integer, and where we define = a and 
b, = b. The integer p determines the representation of 
SU(2); it corresponds to the number of columns in the 
corresponding Young diagram. The spin magnitude is 
simply S = \ p, and the dimension of the representation 
is g = p + 1. 

Now let us add in hole states. The constraint equation 



a f a + 6 f 6 + f\f = p 



(7) 



There are now g = 2p + 1 possible states, corresponding 
to the two classes 



S=\p 

s = \(p-i) 



a ] a + b ] b = p and pf = 

a)a + tfb=p-l and Pf = l 



The simplest such case we shall deal with is p = 2, for 
which g — 5. Explicitly, these states are given by 



1) = ^) » 



l+^) = / t « t l v ) 



1 

V2 

|0) =a t 6 t |V) , 

|-i) = ^(^) 2 |v), |-i> = /t 6 t|v), 

where | is the vacuum for bosons and fermions: a | = 
6|V) = /|V> = 0. 

The 2p + 1 states obeying the constraint of Eq. ([7]) may 
be grouped into a multiplet of the superalgebra OSp(l|2). 
This group has five generators, three of which are the 
familiar SU(2) spin operators: L a — i 6^ <r° b v , with i = 
1, 2, 3. The remaining two generators are non-Hermitian, 
and may be taken to be 



K l = -{x- l fa^ +xfb) 
X 2 = ^(ar- 1 /6t_ a . / t ) 



(8) 



where x is an arbitrary complex number. The relations 
among the generators are 



l€ abc L c I 



(9) 



The algebra of the generators is independent of the pa- 
rameter x. Note that K\ = \a)b = jL + , so K 1 is like the 
"square root" of the angular momentum raising operator 
L + = L 1 + iL 2 . Since \L 3 , KA — \K X , we have that K x 

raises L 3 by half. Similarly, K 2 lowers L 3 by half, and 
functions as the square root of the angular-momentum 
lowering operator L_. 

The Casimir operator is given by 



C = L 2 + e^KK 



(10) 



Acting on the single-site states defined above, C takes 
the value \p(p +1). Generally, one has C = L{L + |), 
where L, which is either integer or half odd integer, is 
the maximum eigenvalue of L 3 , for a given value of C. 
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We call this quantity L the angular momentum. The 
dimension of the representation with angular momentum 
L is g = 4L+1. The addition of angular momenta within 
OSp(l 1 2) is similar to the SU(2) case, except the spacing 
between consecutive L values is \ rather than 1, 



L®U = \L-L'\ 



L-L'\ 



■®{L + L') . (11) 



For example, if the local representation of 0Sp(l|2) on a 
single site is L = 1, then on any link (ij), one can have 



1 3 

l®l=0©-©l©-©2 



(12) 



where the dimensions of the five irreducible representa- 
tions in the product are 1, 3, 5, 7, and 9. The Casimir 
operator for the two-site system is 



c(ij) = 2L{€) ■ + !/•;; f)f. + !/•;, flf, 



(13) 



where 



.4; 



(14a) 
(14b) 



Using the Casimir operator, we can construct projection 
operators onto representations of a desired value of L. 
This can be used to construct a Hamiltonian along the 
lines of AKLT; this program is carried out in Sec IV Al 
The link operator 



V 



= 4b] b 



rftf] 



(15) 



transforms as an 0Sp(l|2) singlet whenever x 2 = —r. 
That is to say [L a , x\j] — [K^ , xj .] = whenever x = 

±ir, where L a = L a {i) and = ^ K^{i) are global 
generators. Thus, if each site is in the L representation, 
there are 2L quanta per site, with p = 2L in Eq.©. 
There are thus AL quanta on each link. In the general 
SUSY AKLT state of Eq.gJ), 2M of these quanta are 
passivated in singlet bonds. Thus, the maximum value of 
for the link is J max = 2L — M, where L = \zM relates 
the value of L, the lattice coordination number z, and 
the integer parameter M. For the L = 1 SVBS chain, for 
example, the wave function is annihilated by projectors 



or Li.i+i — 2 sectors, and 
0, i, and 1. 



onto either of the Z^i+i 
the only remaining possibilities are L i i+1 

However, an inconvenient problem remains. Because 
the generators K ^ are non-Hermitian, this is also the 
case for the projection operators, and, thus, the Hamilto- 
nian as well. Then, we use a "trick" to make a Hermitian 
Hamiltonian from non-Hermitian projection operators as 
performed in Sec|V] In Sec IV Bl we also exhibit a prop- 
erly Hermitian Hamiltonian which has the L = 1 SVBS 
chain with fixed total fermion number an exact nondc- 
generate ground state. Before doing so, we will derive 
the properties of the SVBS chains themselves. 



III. SVBS SPIN CHAINS 

The general SUSY AKLT chain wave function is writ- 
ten as the pair product, 

I*) = IT (4 b li - b l4 + r +rf}fU M |V) • (16) 

i 

This describes a chain in which each site is in the L = Al 
representation of OSp(l|2). The wave function is anni- 
hilated by projectors Pj(i,i + 1) which project onto the 
total link angular momentum J, for M < J < 2M . 

We are interested in computing correlation functions in 
these states. The correlation functions we will compute 
are: 



C sp Jn) = (L(j)-L(j+n)) , 



(17a) 



C ss (n) = ( (a, b j+n - b ■ a j+n ) f] f] +n ) , (17b) 

corresponding to the spin-spin correlation function and a 
"singlet superconductivity" order-parameter correlation 
function. Since our state does not conserve particle num- 
ber, the superconducting order parameter can be non- 
vanishing. Here (O) = | is the normal- 
ized expectation value. A corresponding "triplet super- 
conductivity" correlator, 



C« s (n) 



a j a j+n 
7f \ a 3 bj+n + bj aj +n ) 



V 



bj b j+n 



may also be defined. However, due to the singlet property 
of the SVBS states, we have that C% s (l) = 0. We shall 
compute these correlations on finite chains, which have 
ends, and examine the thermodynamic limit. There are 
some specific properties of edge states in these models, 
in direct correspondence to what is known from AKLT 
chains [H, H^, For example, the edges of the L = 1 
SVBS chain are local L = i degrees of freedom, which 
means that the ground state of a long but finite L = 
1 SVBS chain will exhibit a ninefold quasi-degeneracy, 
with the actual levels arranged into singlet, triplet, and 
quintuplet states, according to^(g)| = 0ffi^ffil. 

Note that the operators whose correlation functions 
are computed must commute with the local constraint 
Expressions such as (/J /| +n ) and 



n a +n b + n f = p. 
(a, bj +n ) vanish identically. 



A. Spin-hole coherent states 

The application of spin-coherent states in elucidating 
the properties of the AKLT VBS states was discussed in 
Ref. [8|. Here we utilize a generalization of the familiar 
SU(2) spin-coherent states, known as spin-hole coherent 
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states Define the state 
\ft,6;p) 

l ) P ® l°)-^ l«) P -i® I 1 ) • ( 19 ) 



3) = — L=(W + w&t - (9/^1 v) 

VP ' 



Here, |n p ) is an SU(2) spin-coherent state with S = ^p, 
and 9 is a Grassmann variable which anticommutes with 
/ and p . The resolution of the identity may be written 
as 



d9 d9e^ eS \h 1 9;p){h,9;p\ =V L= v , (20) 



dh 
-in 



where Pl is the projector onto the angular momentum 
L representation of OSp(l|2). 

Next, consider a general state in the angular momen- 
tum L representation, written as 



VP I 



(21) 



^(at,6t) +V , p _ 1 ( a t ) & t)/t | V ) 



where ip p (a' , ¥) is homogeneous of degree p in and 
w . We then have 



(h,9;p\ip) = ip(u,v,9) 



(22) 



That is, we simply replace a} — > u, — > v, and /t — > 
in the function 0- 



B. Matrix elements 

Now consider the following spin operators: 



m— n— 
fe fe+1 



m— n— 
fe+1 fe 



a a 1 0' 



(23a) 
(23b) 
(23c 



Note that raise (+) or lower (— ) the angular momen- 
tum by AL = i, while T® preserves total spin. Our goal 
is to compute the matrix element 



(i;\f k \ ( t>) = -(V\ij(a > b,f) f£+f+f 



pi 



(24) 



+ P + T fe ' fp] 0(at,6+,/t)|V) 



and to represent it as an integral over spin-hole coherent 
states. We find 

(l>\t k \4>) = J^Jd9jd9e^ s 

xi,{u,v,6)T k {h,6,6)cj>{u,v,6) , (25) 



where 

T k (n,6,6) 



(p+k+l)\ 



p + k + 



j+T^(h)9 (26) 



+ T k -{n)9 + T%{n)< 



MB 



replaces 



% = T° k +f+f + f- /t + f °' //t . (27) 



C. Correlation functions 

With the spin-hole coherent state formalism developed, 
we are now in position to calculate the correlation func- 
tions in the general SVBS chain state. The first step is to 
compute the wave function normalization, which we call 
V (for "denominator"). Using the resolution of unity for 
the spin-hole coherent states, we have 



N 



dp J] 



l n "n+l u n "n+l 



n+l n I ; 



(28) 



71=0 



where the measure is 



N+l 

dp= n 

3=0 



dh- _ 
-r± d9 1 d6, 
in 



x e (2M+l)(8 1 8 1 + ...8 N N ) _ ("29-) 



Note that the site j = and j = N + 1, which are at the 
ends of the chain and have only one neighbor, are in the 
L — \M representation of OSp(l|2) while the bulk sites 
are in the L = M representation. We now expand 



U n V n+l V n U n+1 ' T "n+l "n 



1 n n ' n n+l 



M 



f\ — h ■ h \ M ~ l 
) +Mr( 1 " +1 j {u n v n+1 -v n u n+1 )9 n+1 9 n 



Mr 



1 — h„ • h, 



n "■n+l 



M-l 



+ M 2 \r\ 2 l 1 K ' fln+ < 



( u n v n+l v n u n+l)" n ®n+l 
M-l 

®n "n ®n+l "n+l • 

(30) 



Using 



dh ( 1 — n • h' 

47T 



1 



M + 1 ' 



(31) 



we can now integrate out site j > = 0. The new integrand is 
then the truncated product wave function, starting with 
site j — 1, multiplied by the quantity a 1 + /3 1 9 1 9 1 , where 



6 



a x = 1 and (3 l = M|r| 2 . The form of this expression 
self-replicates. That is, after integrating out sites j = 
through j = n — 1 in succession, we are left with the 
expression a n + (3 n 9 n 9 n . We can now integrate out site 
n to obtain the replication formula, 



X K + P n OJn) X 



, f 2i |2 / 1 n " ' "n+l 



+ M 2 |r| 



1 "n " 

2 

M-l 



A I 



'n w n ™n+l n+1 



2M + 1 1 

m+T a " + m+T Al J + M|r| a ™ 

(32) 



Note that in propagating the expression (a„ + (3 n 8 n 8 n ), 
we may drop the last two terms on the right-hand-side 
(RHS) of Eq.dSDJ). We now have 

v 



' 2M+1 



l n+l 



M+l M+l 

M\r\ 2 



(33) 



When we get to the last site, we have the final result 
C dh A 



V 



Thus, 



4?r J ao N+lJ aa N+l c 

X ( a N+l + Pn+1 ^N+l^N+l) 

(M + l)a N+1 +p N+1 . (34) 
23= (M+l l!P V (.\/>) CiM 



Now we need to compute the numerator for the correla- 
tion function of interest. 



1. Singlet superconductivity correlations 



We define the singlet off-diagonal correlation function 

C SS (») = < &fc+r> - 6 fc a fc+n) /* /*+•» ) ' ( 36 ) 

which is independent of k in the limit of a long chain 
(N — > oo). The operator above, in the language of the 
operators T k studied earlier, is of the form T k=Q on sites 
k and k + n. Invoking Eq. (|26[) . we have 



( a fe kfc+n b k a k+n ) f k f k+n — ► 

(2M + l) 2 ( Ufc Ufe+ „ - v k u k+n ) 9 k 6 k+n 



(37) 



The correlation function may then be written 
C ss (n) = i=(2M + l) 2 .^. 



(38) 



The calculation of the numerator J\f proceeds along the 
same lines as that of 23. Starting with site 0, we generate 
the expression + (3 n 9j9j. When we arrive at site k, 
only the second term on the RHS of Eq. (j3"0|) contributes. 
We then have 



dh, 



«W5W^ (2M+1)V * 



(39) 



x (u, Vj i — v, u, , ) 6, Qj , 



x K ^+1-^^+1 + ^^+11 



1 2M 



M-l 



K v fe+i -« fe « fc +i) 



X K ^/c+n ~ W /c U fe+n) ^fe+l ^fe+n 

/ Mr \ 



When we integrate over site k + 1, we obtain 



M]r| 
M+l 



a /c ( U fc+2 V k+n V k+2 U k+n) u fc+2 "fc+n 



(40) 



We have now replicated the form of the integrand. 
Clearly whenever n is even, the numerator M vanishes. 
For n odd, we obtain 



/M|r| \" 

Pfc+n — \ M + 1 / ™' odd ' fe ' 

The correlation length £(M, r) is then given by 



-l/C(M,r) 



1 / M|r| 



A+ VM + 1 



(41) 



(42) 



M\r\ 



M + | + y / (M+i) 2 + M(M + l)|r| 2 ' 



where 



A, = 



M+| //M+i\ M|r| 



M + l V V M + 1 / M + l 



(43) 



is the largest eigenvalue of the matrix 23 from Eg. ([33 
We can define the s-wave order parameter as 

A = ( K h+i - h «fe+i) fl fl+i ) 

2M(M+ \fr 



(^M(M + l)(l + |r| 2 ) + I + l(M+i)) Z - i(M + l) 2 

(44) 

(A is plotted in Fig. H) 
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0.4 0.6 
sin 2 (77) 

FIG. 2: Order parameter A = ((a n b n +i — b n a„+i) fnft+i) 
in the general SVBS chain. The parameter rj is given by 
r\ — tan -1 \r\. 



2. Spin correlations 

We next turn to the spin-spin correlation function, 
C spin (n) = / L(j) ■ L(j 4- n) ). The spin operator is given 

by L = \h^u illV b vl and is of the type T k=1 . Accordingly, 
Eq. ({2"6")) gives the prescription 



(45) 



Once again, the correlation function is expressed as a 
ratio of TV 'I'D. In the numerator, when we arrive at site 
k, we have the integral 



1 — n, • h 



k "fc+1 



+ M 2 |r| 2 ' 1 



M-l 



7 k"k ^fe+l^fc+l 



( a k+l + P, 



fc+l"fc+l"fc+17 "fc+1 ' 



(46) 



with 



' a k+l\ M(M+h) I W+2 \ ( a k~ 



M + l 



J k+lJ 



,(M-l)|r| 2 / \[3 k 



(47) 



For sites j between k and k + n, we have 



^i[deJdd H e^ M+1 ^ 

47T / 3 I 1 



x (a. 4-/3//,) 

' 1 — n j • h 



1 - • n j+1 



+ M z \r 



2 I |2 



'J "341 



2 

Af-1 



A I 



2 / "i j i+i ;+i 



with 



(48) 



(49) 



where 



/C 



M 



2M+1 



(M+l)(M + 2) l (M _ 1)( Af + 2 )|r| 2 



0; 



(50) 



Finally, we come to site k+n, where we have h k+n -h k+n 
1, and 



X + Afc+rA+rA+n 

1 — fir. , „ ' 72 



1 "fc+n ' 



A/ 



4M 2 |r| 2 



A/-1 



X ®k+n k+n ®k+n+l®k+n+l 



( a k+n+l + Pk+n+l 6 , 



' l~ J k+n+l u k+n+l w k+n+l) ' 



(51) 



with 



vPfc+n+1/ 



= (M + -) 



A/+1 \ / a k+n 







\@k+n/ 



(52) 

For sites / > + n, we propagate by the matrix D from 
the denominator. 

Assuming N — > oo, with n finite but large, we can 
ignore the ends, and we obtain 



C spin (n)=A(Xj\ v ) lnl 

= A (-1)™ e~ |n|/e(M, ' r) 



(53) 



where A is a coefficient, and X K v are the largest magni- 
tude eigenvalues of the matrices K, and D, respectively. 
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0.4 0.6 
sin 2 (?7) 

FIG. 3: Spin correlation length £(M, r) (blue) and supercon- 
ducting correlation length C,(M, r) (red) for the general SUSY 
AKLT chain. 



The spin correlation length is thus given by 



-l/£(Af» _ 2K 
M 



(54) 



(M-l)(M+2)|r| 2 



M 



M(M+l)|r| 2 
(M+A) 2 



In the r - 



limit we recover the result C(n) 



-rr— ,, found for general AKLT chains in Ref. 
[4l| . Note that for r — > oo and M = 1 the correlation 
length vanishes. This is because in this limit the ground 
state is that for the S = | Majumdar- Ghosh model, i.e., 
alternating singlets, for which there are no correlations 
beyond nearest neighbors. For the M > 1 generalizations 
of Majumdar-Ghosh, however, the correlation length is 
finite. The spin correlation length £(Af, r) and super- 
conducting correlation length ((M, r) are both plotted 



in Fig. [3j versus the parameter sin ij 



7(i 



These two correlation lengths coincide at r = (2M +l)/3. 
Especially, when M = 1, they coincide at r = 1. 



IV. RELATION TO QHE STATES 

Here, we discuss analogies between the lowest Landau 
level (LLL) physics and the spin physics in detail. Much 
of our discussion is an extension of the pioneering work 
by Haldane [l(J on the FQHE in a spherical geometry. 

We begin with a discussion about analogies in one- 
particle problem. For The LLL bases are given by the 
monopole harmonics [42j , which form an irreducible rep- 
resentation of SU(2) indexed by the unique Casimir op- 
erator, which is the monopole charge. As is well known, 





LLL physics 


Spin physics 


Space 


External 


Internal 


Quantum number 


Monopole charge 


Spin magnitude 


Basic state 


Hopf spinor 


Spin-coherent state 


Manifold 


Fuzzy sphere 


Bloch sphere 



TABLE I: Correspondences between LLL physics and spin 
physics 



the monopole harmonics in the LLL are constructed as 
symmetric products of Hopf spinors. Mathematically, the 
Hopf spinor is equivalent to a spin coherent state for a 
state in the fundamental (S = |) representation. The 
symmetric products of the spin coherent states give rise 
to higher spin states. In the LLL, the kinetic term is 
quenched and the coordinates of the two-spheres are ef- 
fectively reduced to operators of SU(2) algebra. Such 
manifold with noncommutative coordinates is known as 
the fuzzy sphere and its mathematical structure is equal 
to the Bloch sphere of spin physics. The relations be- 
tween the LLL states and the spin states are summarized 
in Table HI Thus, as for the one-particle problem, there 
are apparent analogies between the LLL physics and the 
spin physics. 



A. Laughlin-Haldane and AKLT states 

Even in many-body level, as briefly mentioned in SecHl 
remarkable resemblances between the Laughlin state and 
the AKLT state have been reported in the work of one 
of the authors @]. On Haldane's two-spheres, particles 
are uniformly distributed to form a rotationally invariant 
incompressible liquid described by the Laughlin-Haldane 
function, 



4h } = 



JV 

m 



ViUjY 



(55) 



where (it, v) indicates the Hopf spinor. Meanwhile, the 
AKLT state is the VBS state made by the SU(2) singlet 
combination of Schwinger bosons [Eq.(|T])]. and, in the 
spin-coherent state representation, is written as 



(M) _ 
AKLT — 



ViUj) 



M 



(56) 



Obvious resemblances may be found between Eqs. (|55[) 
and (|56p . The power m in the Laughlin-Haldane state 
takes even or odd integer depending on the statistics of 
the particles, while M in the AKLT state specifies the 
number of the valence bonds on each site and has nothing 
to do with statistics. Since these two states are "almost" 
mathematically equivalent, their truncated pseudopoten- 
tial Hamiltonians are similarly constructed by the form of 



two-body interactions: the truncated Hamiltonian for the 
AKLT state is given by Eq.Q, while for the Laughlin- 
Haldane state, it is given by 



2.S' 



Kj J' + l 



(57) 



where J* = 25 - m with S = m(N — l)/2. Based on 
the OSp(l|2) supergroup analysis, the SUSY Laughlin- 
Haldane wave function was proposed as 



N 

SLH = Y[( u i v 3 ~ v iUj + r0 l 6 j ) m , 
i<j 



(58) 



where (u,v,6) indicates the SUSY Hopf spinor. In 
Ref . [2(| . r is fixed as —1, but here we take r as a free 
parameter. Extracting the original Laughlin-Haldane 
wave function, the SUSY Laughlin-Haldane state can be 
rewritten as 



(m) 
SLH 



exp I mr 



l<3 



UiVj 



V 1 U 1 



(m) 
LH • 



(59) 



All of the important physics are included in the expo- 
nential factor of Ea.(|59p. and this deformation enables 
us to perform an intuitive interpretation of the SUSY 
Laughlin-Haldane wave function. The denominator of 
the exponential factor, l/{u%Vj — ViUj), represents a p- 
wave bound state of two particles i and j, and the SUSY 
Laughlin state is regarded as a p-wave superfluid on the 
original Laughlin state [28| . By expanding the exponent, 
one may find 



*SLH — LH 



JV 

E 



u^j — ViUj 



■ 4> 



(m) 
LH 



N 
i<j 



UiVj - ViUj 



(m) 
LH 



+ ( m r) N ^l[9 l -A( J] " 

\ . Uj-iVj - Vj-xl 
i \ j:even 



(m) 
LH > 

(60) 

where A in the last term represents antisymmetrization 
over all different choices of breaking particles into pairs, 
and is simply known as the Pfaffian. Hence, the last 
term in Eq. (|6T)|) represents the Pfaffian state proposed 
by Moore and Read HU 



(m) 
MR 



Ui-iVi - Vi-iUi 



(m) 
LH 



Pf 



1 



UiVj — ViUj 



(m) 
LH ' 



(61) 



where all of the particles form p-wave pairings to form 
a bosonic QH state. It is noted that the expression (|60|) 
should be regarded as the expansion about the parameter 
r not m, since the original Laughlin-Haldane function 
itself depends on m. 




-0 




FIG. 4: The single-bond breaking operator annihilates a va- 
lence bond and creates a fermion pair on the nearest-neighbor 
sites. 



B. Physical interpretation of the SVBS state 

Inspired by the similarity between the Laughlin- 
Haldane and the AKLT states, from the SUSY Laughlin- 
Haldane wave function [(Eq. (|58| )]. one may derive the 
SUSY AKLT state, 



(M) 
AKLT 



(«■) 



'i v 3 



ViUj + 



A I 



(62) 



which is the spin- hole coherent representation of Eq.(|3]). 
In the following, we focus on the SVBS spin chain. Just 
as in the SUSY Laughlin-Haldane case, the SVBS spin 
chain state z = 2 is rewritten as 



(M) 
AKLT 



exp 



ViUi+l 



■ $ 



(M) 
AKLT' 



(63) 



where the exponential factor 9i0j / (uiVi+i—ViUi+i) which 
we call the "pair creator" has the following physical inter- 
pretation: it replaces one of the valence bonds between 
sites i and i + 1 by a fermion (hole) pair; this is depicted 
in Fig. [4] The SVBS chain state is expanded as 



(M) 
AKLT 



CM) 
AKLT 



UiV i+ i - ViU i+ i 
2 



■ 4> 



(M) 
AKLT 



UiV i+ i - ViU l+ i 



■ $ 



(M) 
AKLT 



(Mr)^n^ n-n 



U t V i+ i - ViUi + l 



(M) 
AKLT- 



(64) 



We assume here that the total number of sites L in our 
ring is even. The original AKLT state appears as the 
first term in this expansion in powers of the Grassmann 
coordinates. The second term consists of superpositions 
of all AKLT states with one hole pair, the third term of all 
superpositions with two hole pairs, etc. The final term in 
the expansion contains the product 6\ ■ ■ ■ 9l over all sites. 
Its corresponding spin wave function is a superposition 
of two generalized Majumdar- Ghosh states, one in which 
a valence bond has been removed from each even link 
(2n, 2n + 1), and the other where a valence bond has 
been removed from each odd link (2n — 1, 2n). Note that 
each site can accommodate at most one hole (Fig. [S]). As 
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FIG. 5: Graphical representation for the expansion of the 
SUSY AKLT spin state with M = 1 [Eq.flMJ]. At the n th 
term of the expansion, there appears the superposition of 
AKLT states with (n — 1) hole pairs. In particular, the origi- 
nal AKLT state is realized as the first term (r — > 0) and the 
MG dimer states are realized as the last term (r — » oo). 



discussed in fj]] for M = 1 the last term of Eq. (f6"4"|) gives 
precisely the S = \ Majumdar-Ghosh state, 



n-n 



U t V i+ i - ViU i+ i 



(M=l) 
AKLT 



(65) 



where $ A and <E> B correspond to the two dimer states of 
Eq.©. 

Thus, both the Majumdar-Ghosh and Moore-Read 
states appear as the last terms in the expansion of the 
corresponding super wave functions. It is interesting 
to note in this regard that both the Majumdar-Ghosh 
and Moore-Read wave functions vansh when any three 
particles (Moore-Read) or any three neighboring spins 
(Majumdar-Ghosh) coincide, and their truncated pseu- 
dopotential Hamiltonians are constructed by three-body 
interactions [4J|. For the Moore- Read state, 



35 



ff MR = ]t vjPj(i,j,k) 



(66) 



i<j<k J=3(S-m)+2 



with S = h[m(N— 1) — 1], while for the Majumdar-Ghosh 



state 



H MG = Vj5^Pj(i,t + l,i + 2). (67) 



C. More fermion coordinates 

Our construction may be generalized to include addi- 
tional Grassmann coordinates. Introducing two Grass- 
mann species 9i and rji, we write the extended SUSY 



Laughlin-Haldane wave function as 

*slh = li&iVi - v iUj + n BiBj + r 2 rj iVj ) m , (68) 

i<j 

where r\ and r 2 are two free parameters. We may now 
write 



(m) / \ 

SLH = ex P( mr l 2^ 

i<j 



UiVj — ViUj 



(69) 



• exp mr 2 



ViVj 



exp 



i<j 

i<j 



UiVj — VjU 



i u 3 . 



(UiVj - V l Uj) 2 



(to) 
SLH ' 



We have already encountered the first and second expo- 
nents of Eq. ((69|) in the previous analysis, each of which 
represents the p-wave pairing state. The third exponent 
is the newly appeared term, and its exponential factor 
provides (— l) 2 by the interchange of i and j to suggest 
the property of cJ-wave pairing. When we expand the 
third exponent, at the last term, we obtain 



JV 



• exp —mriT2 



E 



nVjViVj 



(UiVj - ViUj)' 



{-mrir 2 ) N/2 S 



(u z Vj - ViUj)' 



(70) 



where Yii^idr/i = Yii d@i Yii ^Vii anc ^ & represents the 
symmctrization operation, which is realized by changing 
all the signs of terms in Pfaffian to be plus, and is known 
to yield the Haffnian, 



S 



1 



(UiVj - ViUj)' 



= Hf 



1 



(UiVj - ViUj) 2 



(71) 



The first and second exponents in Eq. (|69p are expanded 
as in Eq. (f6T)]) to yield the product of two Pfaffians, and 
produce the Haffnian again, 



(mr 1 ) N / 2 {mr 2 ) N ' 2 ¥{ 2 



1 



(UiVj - ViUj)' 



( ri r 2 ) N / 2 m 



(UiVj - ViUj)' 



(72) 



Besides this, there are many cross terms to yield Haffnian 
in the products of expansions of the three exponents. 
Collecting all of the contributions, the last term of the 
expansion [Eq. (|69[) ] is summarized as 

|n***»*S= (m(m-l)™)" /a •*<£>, (73) 



where $hr is the Haffnian state of Haldane-Rezayi [4J 



1 



[UiVj - ViUj) 2 



(m) 
SLH' 



(74) 
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FIG. 6: The operation of the double-bond breaking operator. 
The white circles represent the hole pair of 8i9i+i, while the 
light yellow circles represent the other hole pair of rjirj i+1 . 



which represents <i-wave pairing QH state. We see that 
the Laughlin wave function with two Grassmann species 
[Eq. ([69]) ] is expanded as 



*SLH — LH 



+ {m{m-l)r 1 r 2 ) N,2 Y{6 



(m) 
HR ' 



MR 



(75) 



Intriguingly, with two species of Grassmann coordinates, 
there appear Laughlin, Moore-Read, and Haldane-Rezayi 
states as expansion coefficients. Each of them natu- 
rally appears in the following limits: the Laughlin state 
at ri,r 2 — > 0, the Moore- Read state at r% — > oo or 
r 2 — > oo with r\r 2 fixed, and the Haldane-Rezayi state 
at r\,r 2 — *■ oo. Now, let us move to the discussion of the 
VBS model with two species of Grassmann coordinates. 
The corresponding generalized AKLT state is 




"I" (o^ ^ ^ ^ ^ 

-\m oj vp*^ — — v y — 



+ 



FIG. 7: The graphical representation for the expansion of the 
generalized AKLT spin chain state [Eq, (|77|l ]. The first term 
represents the original M = 2 AKLT state. At the second 
term, the superposition of the AKLT states with one hole pair 
appears. At both third and fourth terms, one may find the 
AKLT states with two hole pairs. At the third term, the two 
holes are generated by the double-bond breaking operation, 
while at the fourth term, they are generated by two successive 
different single-bond breaking operations. At the last terms 
of the expansion, we obtain four states two of which are fully 
dimerized states, and the other two are partially dimerized 
states that are equal to the M — 1 AKLT states. 



(AT) 
AKLT 



JJ(«i«7 - ViUj + ri6i0j + r 2 rnrjj) M , (76) partially dimerized states. An explicit calculation yields 



(ij) 



and, for ID spin chain, it is rewritten 



* AKLT — * AKLT ex P I m r V , 



exp I M r 2 y^^ 



UiV i+ i - ViU i+ i 

ViVi+i 



(77) 



cxp ( — Mr\r 2 



UiV i+ l - 



(u t V i+ l ~ ViUi+lf 



In the following, we concentrate on the case M = 2. 
The factor of the third exponent 9i9i + ir]ir]i + i / (iiiVi+i — 
ViUi+i) 2 is interpreted as the "double-bond breaking op- 
erator" : it annihilates two valence bonds and creates two 
kinds of fermion pairs between i and i + 1 sites [Fig. [6]. 
Then, in Eq. (|77p . there are two types of bond breaking 
operations, one of which is the single-bond breaking op- 
erations performed by first and second exponents, and 
the other is the double-bond breaking operation by the 
third exponent. With this interpretation, we have a nice 
graphical understanding of the expansion of the gener- 
alized AKLT state (see Fig. [7]). As expected from the 
graphical representation, in the last terms of the order of 
(r\r 2 ) L l 2 there appear two fully dimerized states and two 



- AKLT 



(2 n r 2 ) L / 2 n +n (u^ 1 - ^ u *+i) 2 

Vovcn odd/ 

- 2 L+1 (nr 2 ) L / 2 IJ(«<«i+i - «i«<+i) , (78) 

i 

where once again we consider a ring of L sites, with L 
even. Equation ([75)) corresponds to the expression ([75)1 of 
the QHE. The first two terms on the RHS in (J7SJ denote 
the two fully dimerized states, while the last term on 
the RHS represents the two partially dimerized states. 
These fully and partially dimerized states are degener- 
ate zero-energy eigenstates of the three-body truncated 
pseudopotential Hamiltonian, 



H D =J2^2VjPj(i,i+l,i + 2). 



(79) 



,7=2 



The degeneracies may be resolved by adding terms in- 
volving other projection operators to the Hamiltonian 
0. Since the fully dimerized states in Eg. ([75)1 only take 
the spin magnitude J = 1 for groups of three consecutive 
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sites, they are the zero-energy eigenstates of the Hamil- 
tonian, 



H FD =H D +Y, V o P o(h * + 1,* + 2) 

i 



(80) 



while the partially dimerized states are not. 

Comparing the two expressions ([78]) and ([73]) . one no- 
tices the apparent analogies between the fully dimerized 
double-bond states and the HR state. As in the case of 
the dimerized single-bond state and the MR state, they 
share common features such as the truncated pseudopo- 
tcntial Hamiltonians which render them exact ground 
states. For the fully dimerized state, the Hamiltonian 
is given by the three-body interaction form (|80[) . while 
for the HR state, it has a similar form , 

35 

H ™=E E VjPj(i,j,k), (81) 

i<j J=3(5-m)+3 

with S = |(m(iV- 1) - 2). 

The generalization with more fermionic coordinates is 
a straightforward task. With F species of fermionic co- 
ordinates, the SUSY AKLT state is generalized as 

f M 

*AKLT = n(^-^+E r /^/) > ( 82 ) 



and is rewritten as 



(M) _ 
AKLT 



expjM^r/^ 



qfaf 



UiVj — ViUj 



f (<i) 

e X pf-M V r f r' V ^ 

F ei etefet'efet" 

exp( 2M 2^ r/r/zr/"^' 



/</'</" 



(UiVj - ViUj) 3 



exp[(-l) F - 1 (F-l)\Mr 1 r 2 ---r F 



alQla2r\2 



(ij) 



■ $ 



(M) 
AKLT' 



(83) 



As in the previous discussion, we consider the expan- 
sion of the exponentials in Eq. (|83|) . At the first term of 
the expansion, we obtain the original AKLT state with 
S = \zM. The last terms, of order (M F rir 2 ■ ■ ■ r F ) i/2 , 
represent a nearest neighbor RVB state with S = h(zM— 
F). For the SVBS spin chain, the last terms are (fully 



and partially) dimerized states that are degenerate zero- 
energy eigenstates of the three-body interaction Hamil- 
tonian, 
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Hd = J2 E VjPj(i,i + l,i + 2) , (84) 



i J=S+1 



with S = \{2M-F). When M = F, the two degenerate 
fully dimerized states appear in the last terms, and are 
the zero-energy eigenstates of the truncated Hamiltonian 
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with S = \M. 



D. BCS aspects of the SVBS state 

In Sec lIVi we have mainly discussed the property of 
the SVBS state in the two limits r — > , oo and found 
that the M = 1 SVBS spin chain produces the original 
AKLT state at r — > 0, while the MG state at r — ► oo. 
With finite r, the SVBS state contains a finite density 
of hole pairs, and accordingly exhibits superconducting 
properties. This state of affairs is familiar from the BCS 
state, 



BCS > = II 7=1 , 2 i 1 + 9k c{cU) |0> 
V 1 + \9k\ 2 



(86) 



As gk — * 0, the BCS state is reduced to the vacuum, while 
at <7fc — ► oo, it becomes the completely filled Fermi sphere. 
For intermediate gk, the |BCS) describes a state with off- 
diagonal long-ranged order. Then, one may conjecture 
the following correspondences: 



9k <-> r, 



|0) ^ $AKLT, \F) $ 



MG 



(87) 



Interestingly, the BCS state exhibits a duality (S duality, 
in terminology of high-energy theory) with respect to the 
coherence factor, 



9k <-> 1/fffe- 



(88) 



To see this, it is important to notice that the BCS state 
is represented in two ways, 

l BCS > = II -/TTVkT exp {gk c " clk) |0) 



n , 1 



exp 



fcfc 1 ^-*)! ). (89) 



where hk represents the hole operator h) k = c~k, and |0^ 
is the hole vacuum, with hk |0)) = 0, namely, the fully 
occupied Fermi sphere |0^ = |F^. As seen in Eq. 
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the two descriptions in terms of particle and hole opera- 
tors are completely equivalent, and the duality physically 
represents the particle-hole symmetry. The order param- 
eter 



Ai 



9 k 



1 



1 



\9k\ 



9k+9l 



(90) 



manifestly reflects the dual structure of Eq. 1(88]) . The or- 
der parameter thus vanishes in two limits: the weak limit 
gk — > 0, and the strong limit g^ — > oo. It takes its max- 
imum value at the self-dual point \g k \ = 1. The average 
occupancy of the momentum k state, and its fluctuation, 
are given by 



(nk) = 



\9k{< 



l + \9k\ 2 ' 
{{n k -{n k )f) = 



(91a) 



9k 



(i + kl 



(gk + gr'M + gk 1 )' 

(91b) 



The fluctuation, too, is maximalized at the self-dual point 
\gk\ = 1. As the duality is manifest in the BCS state 
and especially between |0) and |F), one may speculate 
a hidden duality between the AKLT state and the MG 
state 



, dual? -T- 
*AKLT < > %G- 



(92) 



Indeed, the parameter-dependent terms in OSp(l|2) 
Casimir operator, Eq. (|13p . are given by 

~ Tr ~ hh " I ~ haj) ti f i> (93) 



which implies a duality 

r <-> 1/r , 



H u 3 



fifj 



(94) 



This is also the case vis-a-vis the truncated pseudopoten- 
tial Hamiltonians for the SVBS states. Physically, this 
duality corresponds to the interchange of VB and fermion 
pair, in which case the SVBS state of Eq.Q is obviously 
invariant under the dual transformation. Though the VB 
and the fermion pair operators possess same antisymmet- 
ric property with interchange of i and j, their squares 
exhibit different properties: the square of the VB is non- 
zero, while the fermion pair vanishes. More typically, we 
cannot naively take the limit r — ► oo in the SVBS state, 
since in that limit, the SVBS state becomes 



KLT 



(ij) 



rjj = 0, 



(95) 



unlike the BCS state. Because of the asymmetric prop- 
erty between VB and fermion pair, the SVBS spin chain 
is not self-dual at the point |r| = 1 and the order param- 
eter [Eq. (f4"4"]) ] takes its maximum value 



| Am „|HV5-2),te^ 



M + l 



(96) 



at 




2M(M + 1) 



(97) 



The expectation values for the boson number rib(i) = 
a\ai + b\bi and the fermion number n,{{i) — fjfi are 
calculated as 



(n b ) = 2M - 1 + 
1 



2M + 1 



V'4M(M + l)(l + |r| 2 ) + l 
2M + 1 



y/4M(M + l)(l + \r\ 2 ) + l 



(98) 



As expected, with increasing |r|, (rib) monotonically de- 
creases, while (rif) monotonically increases. The fluctua- 
tions for the boson number 8n^ = (n^) — (rib) and the 
fermion number Sri{ = (n 2 ) — (rif) 2 are also evaluated as 



Sn h = Sn 2 = x(l — x) 



2M + 1 



^4M(M+ 1)(1 + |r| 2 ) + 1 ' 
and their maximum is Sn^ = 5ni = h at x = or 

r =31 I ' 



(99) 



4M(M + 1) 



(100) 



V. HAMILTONIANS FOR THE SVBS STATE 

In SecslIIII and |lVl we have studied the properties of 
the SVBS state [Eq.(f3])] and its relation to the Abelian 
and non- Abelian fractional quantum Hall wave functions. 
To obtain a better understanding of what physical sys- 
tems the SVBS states describe, we shall in this section 
construct a Hamiltonian for which the SVBS state is a 
unique ground state. 



A. Generic truncated pseudopotential Hamiltonian 

As mentioned in Sec|TTl the SVBS state [Eq.(gJ] is 
invariant under OSp(l|2) transformations generated by 
the parameter-dependent generators, L a and K^, when 
x 2 = —r. Taking advantage of this symmetry, it is pos- 
sible to construct pseudopotential Hamiltonians for the 
SVBS states with arbitrary values of the parameter r. 
Truncated pseudopotential Hamiltonians for the SVBS 
states [Eq.Q] are constructed by following the similar 
methods of the original AKLT model. The superspin op- 
erator on site i, Li — i(a\ai + b\bi + fj /,), acts the SVBS 
state to yield the eigenvalue L = ^zM. The z component 
of the bond superspin Jf- = L z (i) + L z (j) = \{a\ai + 
Ojaj — b\bi — bjbj) counts the difference between the pow- 
ers of a and b in the SVBS state [Eq. Q] , and the maximal 
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value of J z reads as J^ ax = (z - 1)M = 2L - M. Since 
the SVBS state is invariant under the 0Sp(l|2) trans- 
formation, the maximal magnitude of bond superspin 
is equal to that of its z-component, i.e., J max = J^ ax . 
Thus, the SVBS state does not contain any 0Sp(l|2) 
angular-momentum components larger than J max and is 
a zero-energy ground state of the truncated pseudopo- 
tential Hamiltonian, 



2L 



H 



E E ^ F ./fc'), 



(101) 



(ij) ,/=J max +i 



where Vj are positive coefficients. Pj(ij) is the projec- 
tion operator made by OSp(l|2) Casimir operators, 



2L 

n 



J{J +\)- J'{J' + ±) 

2K A {i)K A {j) + 2L(L + |) - J' (J' + i) 



J(J + i)- J'(J' + 5) 



(102) 



which projects to the two-site subspace of the bond su- 
perspin J. Here, we have used K A (i) — K A (j) = L(L + 
i) with K\ — L'l + e^K^Ky. Apparently, the projec- 
tion operator [Eq. (|102p ] is OSp(l|2) invariant, and hence 
the truncated pseudopotential Hamiltonian (jlOip as well. 
Following similar discussions in the AKLT model, one 
may prove that the SVBS state is the unique zero-energy 
eigenstate of the Hamiltonian (jlOip . 

As an explicit example, it would be worthwhile to 
demonstrate the truncated pseudopotential Hamiltonian 
for the L = 1 SVBS spin chain. With the OSp(l|2) de- 
composition rule (HU), Eq. (|10ip becomes 

ffchain = E( y | P | (*' 1 + l ) + ^ Mh * + 1)) 

i 

= E (H^ - 7Vs)(K A m A ( l + 1)) 4 

i ^ 

+ -(y 2 -w i ){K A {i)K A {i + i)Y 

+ l( m - 7Vi )(K A (i)K A (i+l)) 2 
+ ^(5V 2 + 63Vs) K A (i)K A (i + 1) + Va 



(103) 



In the special case V 2 = 7V3/2, the first term on the last 
RHS in Eq. (|103p vanishes, and (|103[) is reduced to 



with bond superspin | or 2, 
F im (i,i + 1) 



J=0, 1/2,1 



8(K A {i)K A {i + l)) 3 + 28(K A (i)K A (i + l)) 2 

(105) 



fi3 45 
+ -K A (i)K A (i + l) + - 



However, Hamiltonian (jlOip cannot correspond to that 
of any physical system, since it is non-Hermitian 1 be- 
cause of the term t^ v K^K v , as mentioned in Sec|TTJ To 
obtain a physical Hamiltonian for which the SVBS state 
is its unique ground state, one can replace the Hamilto- 
nian (jlOip by the following form: 



2L 



F = E E VjF*j(ij)Fj(ij) 



(106) 



(ij) .!=.]„ 



in which Vj > just like in Eq. (|10ip . Here we would 
like to make several comments on some properties of 
the Hermitian Hamiltonian. First, the definition (|106p 
is a natural generalization of the original pseudopoten- 
tial Hamiltonian, since, if the projection operators were 
Hermitian, with the property Pj = Pj, Eq. (|106p would 
be reduced to the original form (jlOip . Second, unlike 
the non-Hermitian Hamiltonian (jlOip . Eq. (|106p is not 
OSp(l|2) SUSY invariant, because the Hermitian con- 
jugate of the OSp(l|2) Casimir operator contained in 
is no longer invariant under the original OSp(l|2) 
transformation. Consequently, the excitation spectrum 
of the Hermitian Hamiltonian is not SUSY invariant, 
even though the ground state remains is a SUSY sin- 
glet. Third, Hamiltonian (|106p does not preserve the 
total fermion number Nf = fjfi since the Casimir op- 
erator (K A (i) + K A (j)) 2 contains pair-creation terms of 
fcrmions, as shown in Appendix [Bj This is in agreement 
with the fermion number fluctuation in the SVBS state 
[Eq. ([!])]. Physically, such a pseudo-potential Hamilto- 
nian describes some interacting electron system coupled 
with a superconducting bath, which provide a particle 
bath through proximity effect. 

Since Pj(ij) Vj(ij) is always non-negative, it is 
straightforward to prove that H = for a state |G) if 
and only if ¥j(ij) = for all sites and all J max < J < 2L. 
Consequently, if the SVBS state is the only zero-energy 
eigenstate of Hamiltonian (jlOip . it must also be the 
unique ground state of the Hermitian Hamiltonian (|106D . 
One can then prove the SVBS state to be the unique 
ground state of Hamiltonian (|106p following exactly the 



ffchain^E F f P 3/2©2(M + l), (104) 

i 

where P3/2©2 is the projection operator onto the space 



[1] Though the Hamiltonian H101H is non-Hermitian, its eigenvalues 
are still real. Recently, the study of such non-Hermitian Hamil- 
tonians with real eigenvalues has attracted much attentions 46] , 
and the present Hamiltonian would be an interesting example. 
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same procedure as AKLT's original work [|| Q. We will 
leave the detail of this proof as the task of Appendix O 
Here, we sketch the proof for L = 1 superspin chain. Let 
Y&g) De a ground state of Hamiltonian (|106p and satisfy 
the equation 



; > = o. 



Then, 



(* G |-ff|*G>=0 
=> (yG\^j(ij)Fj(ij)\it> G ) J>Jv 
^P./(*j)|* G ) J>Jmax =0, 



= 



(107) 



(108a) 
(108b) 



where in the second arrow [Eq. (|108ap ] we have used that 
Vj in Eq. frOTI) satisfy Vj > 0. Meanwhile, if |* G ) is 
annihilated by the projection operator, i.e., if 



Fj(ij)|*G> 



J>J„ 



0, 



(109) 



then it immediately follows that H\^g) — 0. Thus, the 
condition (|109|) is the necessary and sufficient condition 
such that the | ^g) is the ground state of the Hamiltonian 
(|106l) . We use the condition (|109[) to show |* G ) is the 
unique ground state of the Hamiltonian. For L = 1, the 
condition (|109|) is given by 



•»(»,*+ 1)|*g) =P 2 (m + 1)|* g ) =0. 



(110) 



As we assumed, there is superspin 1 on each site of the 
chain, and therefore, if the two superspins on sites i and 
z + l did not combine a OSp(l|2) singlet, their bond su- 
perspin inevitably would exceed J ma x = 1 due to the 
OSp(l|2) decomposition rule [Eq. ([T2"|) ]. This observation 
holds for bond superspins on arbitrary two neighboring 
sites. Then, on any two neighboring sites, the bond su- 
perspin should form a OSp(l|2) singlet, and the "bulk" 
ground state is given by the products of neighboring 
OSp(l|2) singlet states. Hence, with periodic boundary, 
it is apparent that the SVBS chain state [Eq. (Til)]) ] is the 
unique ground state. With open boundaries, there are 
ninefold quasi-dcgcncratc ground states corresponding to 
directions of the superspins on two ends, 

L-l 

I*g>^ = vio • n - & H+i + rftfU) ■ 4, L \o), 

(111) 

where \i,v = a, 6, /. These ninefold quasi-degenerate 
states generally take different expectation values for local 
observable A, 



(* G |A|^ G ) 



(112) 



However, as in the original AKLT case Q, the differ- 
ent energy eigenvalues converge in the infinite chain limit 
as we shall see below. Suppose the length of the chain 



N (from site to site N), and A takes its support in 
{I, . . . , N — 1} (I <C N). First, we discuss the integration 
of the numerator of Eq. (| 1 1 2(1 from one end (site 0) to site 
I. The inner products of the superspin states at site are 
denoted as 



(113) 



The self- inner products of ito, i>o, and 9q correspond to 
(ao,A),7o) = (i>f>0)> (§;-■§! 0); and (0,0,1), respec- 
tively. The integration from site j to site j + 1 induces 
the transformation: 



(114) 



The three eigenvalues of the transfer matrix are given by 
A± = (3± \/9 + 8|r| 2 )/4 and -|, and then, at I oo, 
the product of the transfer matrices provides 





( "i+l 













(115) 



Then, if there is u or v at site 0, we have a factor 
(l-0/0f) A + fl / 2 (A+ -A_) at site I, while if 9 , we have a 
different value (1 — OiO'i) A + A_/(A+ — A_), but the results 
only differ by the scaling factor, and such difference is not 
relevant to (A) since the scaling factor is canceled be- 
tween the numerator and the denominator in Eq. (|112p . 
Thus, the integration is not relevant to directions of the 
superspin at site in the infinite limit. The integration 
from the other end (site N) to site N — I gives same con- 
sequence. Then, regardless of directions of superspins on 
boundaries, the expectation value of any local observable 
provides a unique value 



(116) 



and, in this sense, the ninefold quasi-degenerate SUSY 
ground states converge to the unique ground state on 
infinite chain. 



B. Another Hamiltonian for fixed total fermion 
number 

In this subsection, we will show an alternative Hamil- 
tonian for the simplest L = 1 case, which is not con- 
structed from the OSp(l|2) Casimir operators but has 
the advantage of respecting fermion number conserva- 
tion. Motivated by the three-site Hamiltonian known 
for Majumdar- Ghosh spin chain [21], here we construct 
a Hamiltonian with both two-site and three-site terms, 
for which the projection of the SVBS state [Eq.([4])] to 
a fixed total fermion number is a unique ground state. 
Such AKLT states with fixed fermion number have ap- 
peared in each order of the expansion of the SVBS state 
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FIG. 8: (a) Schematic picture of a spin configuration of the 
SVBS chain. The blue or orange sites are spin-1 and spin- 7; , 
respectively. k n and l n label the last site of each S = \ (S = 
1) segment, (b) Schematic picture of the ground states of 
Hv in Eq. (|117|l . Each solid line stands for a nearest-neighbor 
singlet pair. The spin of S = 5 sites are free except for the 
neighbor sites of the S = 1 segments, (c) Schematic picture 
of the ground states of Hv+Hu- The S = 1 sites form AKLT 
state (terminated by a S = \ site) and the S = \ sites form 
dimerized MG state. For a fixed configuration of S — 1 and 
S — i sites, the ground state is unique, (d) The effect of the 
hopping term H t , which hops a nearest-neighbor singlet from 
i,i + 1 link to i + 1, i + 2 link, or vice versa. 

as seen in Sec lIV Bl For simplicity, we will focus on the 
M — 1 case, i.e., a chain with S — 1 or S = \ on each 
site. We will first write down the form of the Hamiltonian 
before analyzing the physical meaning of each term. 

i 

H v = ^(F|P|(i,i+l) + F 2 P 2 (i,i + l)) , 

i 

H v = Y^U>P i (i,i+l,i + 2) , 

i 

H t = -tJ2(^i,i+i^l+i,i+2 + h - c -) ( 117 ) 

i 

with Pj(i, i + and Pj(i, i + 1, i + 2) the two-site and 
three-site projections to total SU(2) spin J states, respec- 
tively, and A M+ i = f-f- +1 (a,-6,+i - Mi+i) the annihi- 
lation operator of a Cooper pair. It should be noticed 
that the Hamiltonian is defined in the Hilbert space sat- 
isfying the constraint a\ai + b\bi + fj fi = 2, Vi. The coef- 
ficients V2, Ua, t are all positive. The chemical-potential 

term — /x^j fjfi determines the fermion number in the 
ground states. 

To understand the ground-state property of Hamilto- 
nian pi7[) . we start from the interaction terms Hy + Hjj . 
Since Hy + Hjj preserves the fermion number = fjfi 
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on each lattice site, one can focus on studying its matrix 
element within a subspace defined by fixed eigenvalue of 
n 1 -. For any given configuration {n^ 1 } , the ID chain 
can be viewed as consecutive staggered sectors of spin-1 
and spin-^ chains, as shown in Fig. [8^a). When {nj 1 } 
satisfies 




for k„ < i < l n 

1 for l n <i < k n +i , 



with n € {1, ...,M}, the chain consists of M spin-1 
chains with lengths /„ — k n and M spin-i chains with 
lengths k n +i - l n - (Here k M +i = h.) 

Now we consider the effect of Hy and Hjj on such a 
spin chain. Firstly, the two-site projector P2(i,i + 1) is 
nontrivial only when there are no fermion on the two 
sites (i, i + 1), because the total spin is automatically 
smaller than 2 if there are one or two holes on these 
two sites. Therefore, the V2 term in Hy is an AKLT 
Hamiltonian acting on the disconnected spin-1 segments 
k n < i < l n . Thus we immediately know that the V2 
term takes the minimal eigenvalue of zero if the spin-1 
segments k n < i < l n a- re a U spin-1 AKLT spin chains. 

Second, the two-site projector ^3/2(1, i+1) is nontrivial 
only when there is one fermion on the two sites + 1), 
i.e., vlI + n^ +1 — 1. For these sites, the requirement 
-P3/2 + 1) =0 leads to singlet pair between the free 
S = I spin at the end of the AKLT spin-1 chain and 
the neighbor spin-i site. This requirement automatically 
fixes the length of each spin-^ segment k n +i — l n to be 
> 2. Other spin-i sites which are not neighbor of spin- 
1 site are not affected by Hy. In summary, the spin 
configuration with vanishing eigenvalue of Hy is shown 
in Fig. E](b). 

Third, the three-site projector £3/2(2, i+1, i+2) is non- 
trivial only when there are one or three fermions on the 
three sites (i, i + 1, i + 2). When there are one fermion 
on the three sites, it can be proved that any spin con- 
figuration which satisfy Hy = also satisfy Hjj = 
automatically. Thus we only need to consider the effect 
of Hjj on the sites with three fermions, i.e., three consec- 
utive sites with = nf +1 = n'- +2 = 1. In other words, 
Hjj is exactly the Majumdar-Ghosh Hamiltonian for the 
S = i segments. As known from the work of Majumdar 
and Ghosh, the ground-state requirement Hjj — can 
only be satisfied by the two valence bond solid states, 
with spin singlet pairs between each two nearest-neighbor 
sites. Moreover, the connect condition to the S = 1 seg- 
ments will pick one of the two VBS states, as shown in 
Fig. E](c). (Also, the length of each S = \ segment is au- 
tomatically required to be even, in order to form singlet 
pairs.) 

In summary, the ground state of interaction terms 
Hjj + Hy is uniquely determined for a given distri- 
bution of S = 1 and S — \ sites. Now we con- 
sider the effect of the hopping term H t . The oper- 
ator A^i + i annihilates a singlet pair and creates two 
fermions on i and i + 1 sites. Thus Aj^+iAt, 1 i+2 flips 
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a singlet from i,i + 1 link to i + l,i + 2 link. No- 
tice that A] +M+2 = (al +1 fol +2 - & l+i«l +2 ) fi+zfi+U w e 
know that the term A^+iA]^ i+9 has nonzero matrix 
element only if n^ +1 = n^ +2 = 1, nj 1 = 0. In other words, 
Ht only acts on the interface sites between S — 1 and 
5 = ^ segments. Moreover, in the ground-state mani- 
fold of Hy + Hjj, the effect of H t is simply hopping of a 
nearest- neighbor singlet, as shown in Fig. [8] (d). From 
this picture we know that Ht preserves a ground state of 
Hy + Hjj in the ground-state manifold. Consequently, 
H t lifts the degeneracy of the ground-state manifold of 
Hy + Hij. The lowest energy state determined by H t 
for a fixed total fermion number is obviously the equal 
weight superposition of all the spin configurations satis- 
fying Hy + Hjj = 0, which is exactly the SVBS state 
[Eq.((4|)] projected to a fixed fermion number, 

\G N ) =PnY[ {a\bt +1 - b\a\ +1 + /t/t)| > (118) 

i 

It should be noticed that \Gn) is nonvanishing only when 
N is even, otherwise the ground state cannot be a spin 
singlet. As the last step, the fermion number N for 
which the state |Gjv) has lowest energy can be tuned by 
the chemical-potential term — ft^t f} fi- It is possible 
that for some fj, the ground state contains odd number of 
fermions, which thus cannot be SVBS state. 

In conclusion, we have shown that Hamiltonian (|117p 
has the SVBS state [Eq. (|118p ] as its unique ground state, 
as long as t, V3/2, V2, C/3/2 > and the chemical poten- 
tial is chosen properly so that the ground state has even 
number of fermions. We have also confirmed this fact 
numerically by diagonalizing Hamiltonian (|117p for up 
to five sites with periodic boundary condition and cal- 
culating the overlap between the numerical ground-state 
wave function and the projected SVBS state [Eq. (|118[ l]. 
Within numerical accuracy, the ground state of Hamilto- 
nian (|117p for even total fermion number is unique and 
always given by the SVBS state [Ea.([TTg])]. 

VI. CONCLUSIONS 

In conclusion we have constructed the supersymmet- 
ric generalization of the valence bond solid states. In 
one dimension, these SVBS states smoothly interpo- 
lates between the integer and half-integer VBS states, 
and they represent superconducting valence bond liquid 
states. We also constructed microscopic Hamiltonians 
for which these states are the exact quantum ground 
states. We show that the SVBS states are analogous 
to bosonic Pfaffian states of the quantum Hall effect, in 
precisely the same sense as the analogy between the VBS 
states and the Laughlin quantum Hall states. Our work 
also provides a precise mathematical realization of some 
ideas in strongly correlated systems, in the sense that 
the doped valence bond liquid states are naturally super- 
conducting, and that the superconducting states can be 



obtained from a symmetry rotation, in our case a super- 
symmetric rotation, of the quantum antiferromagnetic 
ground states. For the future, we propose to focus on the 
two- and higher-dimensional versions of the SVBS states. 
Given the analogies between the SVBS states and the 
Pfaffian states in the quantum Hall effect, it would also 
be interesting to explore the possibility of non-Abelian 
statistics of the elementary excitations. 
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APPENDIX A: OSp(l|2) AND SU(2|1) ALGEBRAS 

Here, we review basic properties of OSp(l|2) and 
SU(2|1) algebras with emphasis on their relation to 
Schwinger boson and slave fermion formalism. The 
OSp(l|2) algebra consists of five generators, La = 
L a (a = 1, 2, 3) and L^/j, = 9\, O2) that satisfy 

[Lai Lb] — i e abdL c , 
[LaiLp] = — (<J a ) uti L Vl 

{L^,L V } — ^(e<7 a )^L a , (Al) 

where a a are Pauli matrices, and e is the 2x2 antisym- 
metric matrix e = 102 . Equation (|A1[) suggests that L a 
transform as SU(2) vector and SU(2) spinor. The 
Casimir operator for the OSp(l|2) group is given by 

C = L A L A = L a L a + e^L^L v , (A2) 

and its eigenvalue is L(L+ i) with integer of half- integer 
L. L is referred to as superspin and characterizes the irre- 
ducible representations of OSp(l|2). The dimension of ir- 
reducible representation with superspin L is 4L+1, 2L+1 
of which is the SU(2) spin L representation, and the re- 
maining 2L is SU(2) spin L-|. Specifically, the OSp(l|2) 
fundamental representation L — h is three-component 
spinor and the corresponding OSp(l|2) generators are the 
following 3x3 matrices: 

l a = -h °V l» = -( ° , A (A3) 
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with ri = (1,0)* and r 2 = (0,1)*. The irreducible de- 
composition for superspin representations is given by 



and the superspin magnitude corresponds to the half of 
the total particle number 



L®L' = \L-L'\®\L-L'\ + 1/2®\L-L'\ + 1®---®L+L'. 

(A4) 

Unlike the SU(2) decomposition rule, the superspins on 
the RHS differ by \. 

The SU(2|1) or OSp(2|2) algebra consists of eight gen- 
erators; L a , [OSp(l |2) generators], and T that sat- 
isfy 

[L a , Dfj] = -{<Ja)u^D Vl 

{D^D U } = --(e<7 a ) M „L a , 

{Lfj,, Du} = --e^F, 

[L a ,T] =0, 
[L lt) T] = -D lt> 
[D„T] = -Lr 



(A5) 



With Eq. (|A3[) . the simplest matrix realization for 
Eq. (|A5|) is given by 








1 = 




(A6) 



As the Schwinger particle used in the SU(2) spin for- 
malism, slave fermion is introduced in the superspin for- 
malism. We denote Schwinger bosons as SU(2) spinor 
6 M = (a, b), and slave fermion as SU(2) singlet /, and 
they satisfy the commutation relations: [a, af] = [6, tf] = 
{f,P} = 1- The SU(2|1) operators are represented as 

L, = ^l^ = ±(blf + e^pb„), 
D fl = ^d^ = ±(-blf + e„ v pK), 



T = -tfyil) = a) a + tfb + 2pf, 



(A7) 



and are not 



where tp = (a, 6,/)* = (61,62,/)* 

Hermitian in the conventional sense, while with the def- 
inition of the superstar conjugation J 



(/*)* = "/, (/1/2)* =flfl, 
they become pseudo-Hermitian operators 



4 



(A8) 
(A9) 



L=i(at a + &t& + /V). 



(All) 



We introduce a complex parameter x to define one- 
parameter family of fermionic generators made by L M 
and D.. 



K lt = -^(L lt -D ll ) + ^(L lt + D lt ) 



(A12) 



At x = 1, Kn is reduced to L„, and at x = i, = iD^. 
Though K/j, depends on the parameter x, interestingly, 
L a and satisfy the parameter independent OSp(l|2) 
algebraic relations 

[Lai Lb] — i^abcLc, 



[La,,K(i] = -^{ a a)utiK v , 

{Kfj,,K v } = -(eo-ajftvLa. 
The Casimir operator is given by 



(A13) 



= Lt 





1 \ 2 






(§- 


2x ) e ^ v ^ u + 


(5- 


2x J 



(A14) 

which, in the slave fermion representation, is expressed 
as 



Ki = 



ta + 6t& + /t// o t + &t& + /t/ 1 



(A15) 



Again, the parameter x does not appear in Eq. (|A15P 
and the eigenvalues of the Casimir operator are given 
by L(L + 1/2) for any of the one-parameter family. 



APPENDIX B: SUSY SPIN-SPIN 
INTERACTIONS 

Here, we discuss several properties of the OSp(l|2) 
spin-spin interaction 



(The detail definition of the superstar conjugation can be 
referred to Ref. (47].) In the slave fermion representation, 
the OSp(l|2) Casimir operator (|A2[) is rephrased as 

C^ ta + bt 2 b + /f/ ( ata + bt 2 b + ^ + l), (A10) 



K A (i)K A (j) = L a (i)L a (j) + e^K^K^j), (Bl) 

where i and j represent the sites on which superspins are 
defined. Since the SUSY spin-spin interaction operator 
commutes with the superspin-magnitude operator Li — 
1/2 (at a + b'b + f j 'f)i, the SUSY spin-spin interaction 
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does not change the magnitude of the superspin on each 
site. The bosonic spin-spin interaction part of Eq lBl|) 
gives the SU(2) spin-spin interaction 

L a (i)L a (j) 

(B2) 

while the fermionic spin-spin interaction part of Eq. (|B1|) 
provides 

<^,iuK ,,{1) K v {j) 

= (f + h) e ^ L ^)L»(j) + (l - ^) ^ D ^) D u{3) 

+ (l + h) (f " ^) e ^0M*)A^) + D»(i)L v (j)), 

(B3) 

and, in the Slave fermion representation, expressed as 

+ \{a\a 3 + b\bi)f)fi + \{a) ai + b)b % )f}f 3 . (B4) 

The first two terms on the RHS in Eq. (|B4|) are particu- 
lar interactions existing in the 0Sp(l|2) spin-spin term. 
They violate the total fermion number conservation, and 
represents hole-pair annihilating and creating interac- 
tions. Besides, they are not Hermitian even at x = 1. 
[However, at x = 1, they are pseudo-Hermitian with the 
definition of the superstar conjugation (|A8|) .] The last 
two terms represent interchange of fermion and boson 
between i and j sites. Meanwhile, the SU (2 1 1) spin-spin 
interaction is given by 

L a (i)L a (j) + e^L^L^j) - e^D^D^j) - ^T(i)T(j) 
= ^a\ajb^b t + iaja^jfr., + ^(a f a - 6 t &) l (a t a - & f % 

+ i(a].a, + b^fUi + \(4*j + 4h-)f}fi 

- i(a+a + tfb + 2/t/) l ( a t a + b% + 2pf)j, (B5) 

and is the component of the SUSY t — J model Hamil- 
tonian. It should be noted that the particular hole-pair 
creating and annihilating terms in Eq. (|B4[) do not exist 
in the S U (2 1 1 ) spin-spin interaction. 

Though the OSp(l|2)-invariant spin-spin interaction is 
not Hermitian, its eigenvalues are real and do not de- 
pend on the parameter x. Indeed, with two-body oper- 
ator KA{i,j) = K A {i) + K A {j), the 0Sp(l|2) spin-spin 
interaction (|B1|) is simply rewritten as 

K A (i)K A {j) = \K A {i,j) 2 -\K A {if-\K A {jf , (B6) 



and its eigenvalues are 

E = \j{J + \) - + 1) - ^(Lj + \) , (B7) 

where J, Li, and Lj are the Casimir indexes for K A (i,j), 
K A {i), and K A (j), respectively. 

One may confirm above features with a low energy 
example. The two-body states |j, J3) made by Li = i 
and Lj = i, carry the OSp(l|2) Casimir indexes J = 
0^, 1 by the decomposition rule (|A4j) . The J = sector 
consists of 

\Q 1 0) = (a\b]~b\a]-x 2 f}f})\Q) . (B8) 

This is the 0Sp(l|2) singlet state, and is the "component" 
of the SVBS state [Eq.®]. The J = \ sector consists of 

l^) = (4/j-//«])|o), 

\\,Q) = {a\b]-b\a)-2x 2 f}.f})\Q), 

\\-\) = Q>\f)-fWM- (B9) 
Similarly, the J = 1 sector consists of 
|l,l>=atoj|0>, 

|i^) = (4/j + /H)l°)' 

|1 ; 0> = (a|&t + 6| a t)|0> 5 

\h-\) = (p\f] + ftb))\Q), 

|l,-l) = 6j&t| >. (BIO) 

Equation (|B7j) suggests that J = 0, J = \, and J = 1 
sectors carry eigenvalues E — — hi~ j, an d j, respec- 
tively. By applying the OSp(l|2)-invariant spin-spin in- 
teraction operator to these states, one may confirm such 
parameter-independent eigenvalues are obtained. The 
parameter dependence appears only in the eigenstates 
|0,0) and |±,0), as found in Eqs.jBl} and fB9|) . 

APPENDIX C: PROOF OF THE SVBS STATE AS 
UNIQUE GROUND S TATE OF HAMILTONIAN 
((1061) 

In this appendix we will prove that the SVBS state 
[Eq.Q] is unique ground state of Hamiltonian (|106p . The 
procedure of this proof is a straightforward supersym- 
metric generalization of AKLT's original work [a, 0|- To 
finish the proof, we need to consider the open-boundary 
condition. The boson and fermion can be written in a 
OSp(l|2) spinor tpi = (a,;, bi, fi), and the SVBS state can 
be written as 

|SVBS) = J] {4^ v 4 + i,u) M |0>, (CI) 
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where 



' 1 s 
-10 
^ r, 



(C2) 



For a open-boundary chain with length L, the defini- 
tion needs to be modified by 



/ M \ L-l 

|SVBS;{/^ t }) = ni n(i r 4,r 




1=1 



M 



|0> = n^ t |0) (C3) 



in which fl,, 



"fi B ut — *"^i^i---^m\v\vi---v m ^ symmetric un 
der the permutations (pi^2 ■ ■ • Mm) and (^1^2 ■ • ■ ^m)- 
In other words, the state |SVBS; {^i s , v t }) carries the 
OSp(l|2) representation 4j-<8> 4jr- For the open-boundary 
system, we have the following lemma: 

• Lemma 1. On an open-boundary chain with 
length L, if a state |*) satisfies = 

0, V i = 1, 2, L - 1, iV = M + |, . . .', 2M, then 
the state is a superposition of the SVBS states (|C3[) . 

1. e., 



(C4) 



Lemma 1 can be proved by induction as follows: 
(1) In the two-site case L — 2, the states in the Hilbert 
space are classified by the superspin as 



1 

M®M = oe-ei 



2M. 



(C5) 



,2M 



The requirement P$|#) = 0, V JV = M + ±, 
requires the state to stay in the sub-Hilbert space of 
0®|ele---9M which has a dimension of 1 + 3 + 
. . . + (4M + 1) = (2M + l) 2 . On the other hand, the 
(2M + l) 2 states |SVBS, {fi s , v t }) are linearly indepen- 
dent and satisfy the constraint. Consequently, the states 
|SVBS, {fi s , ft}) span a complete basis of the ground 
state Hilbert space. In other words, the lemma 1 for 
L = 2 is proved. 

(2) An arbitrary state in the Hilbert space 

of a length L + 1 chain can always be expanded as 



|*> lli+ l = £„, TO |»>i L *r 



1) with |n) and 



|?7i) L+1 an arbitrary set of basis states for the Hilbert 
subspace of the first L sites and that of the last site. 
By an SVD decomposition of the matrix ^ n m, one can 
always obtain the form 



l*W = £ A *l^>i,^l^ + i 



(C6) 



Hilbert space of the L + 1th site. The coefficients Xk > 0. 
If P^+i l+1 = for i = 1, 2, . . . , L — 1, we have 



= norm 



^l(Wk\^ i+1 F^ i+1 \W k ) 1 



E i2 

k 



=> 0=\W k \ L = A^\SVBS,{»s,»t}) 1 , L (C7) 

The last step is inductive, assuming the result holds true 
for a system of L sites. Thus the state 1 is written 

as 



L+l 



^^r^r^^lo) (C8) 



in which = Uti 4+i,«r» IL=i |°>- Thc in 



dices (f 5 ,(Jfc) carry the representation 4^ ® 4^, which 
can be decomposed into irreducible representations as 
M(g^ = o©i©---©M. Such a decomposition can be 
expressed as 



M 



N 



= E E *i 



0) (C9) 



7V = 0n=-iV 



in which CX fc , " = ~N,-N + 3,..., N are the 3j- 
symbols carrying the representation of 41- <8 4p (8 A. Thus 
in the state C^fi^ t O^+ T ||0), the sites L and L + l 
carry the representation 4j- ® N (& 4J-. Thus we know 
that the maximal total OSp(l|2) "angular momentum" 
of these two sites is M + N. Consequently, the require- 
ment l+1 = 0, iV > M requires that only 
N = n = terms are nonzero in Eq. (|C9[) . In other 
words, the state can be written as 

= Ko Tl cX k nl±n L a + T ]\o) . (cio) 



Moreover, the coefficient C^* 1 can be expressed as 



no 
a 1 



°00 



S 



.8=1 



(Cll) 



in which S [■ ■ ■] stands for symmetrization over the in- 
dices {vt} and {cfc}, respectively. By the definition 
of fi^, t in Eq. (|C3p . it is straightforward to show that 

C^'^uMtn = Kfn 1 , so that 



I*)i,l + i = FoT'^+'lO) = F^|SVBS,K,ti}) 



(C12) 



where |W r fc) 1 L are orthogonal states in the Hilbert space 
of a length-L chain, and \Sk) are orthogonal states in the 



In summary, we have proved lemma 1 by induction. By 
making use of lemma 1 , it is straightforward to prove that 
the SVBS state [Eq.(5J|] to be the unique ground state of 
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Hamiltonian (| 106|) . First of all, it is easy to see that for 
any physical state = EiE^M+i^v • 

norm(P^ +1 |*» > 0. Since the SVBS state [Eq.Q] 
satisfies J?|SVBS ) = 0, we know that it is a ground state 
of Hamiltonian (| 106|) . On the other hand, if there is 
another state |G) satisfying H\G) = 0, we have 

norm(Pf 4+1 |G)) = 

p 5+i| G ) = 0. y h VM < N < 2M.(C13) 

Consider a chain with L sites and periodic boundary con- 
dition. According to lemma 1, the conditions |G) = 
for i = 1,2, 1 lead to 

\G) = A^n^ t \o). 

In the same way as has been used in the proof of lemma 



1, the coefficient A^' Vt can be decomposed into different 
irreducible representations as 

M N 

A ^ t = J2 FNnCft? • (C14) 

N=0n=-N 

Applying the condition x I G) = to the state | G) = 

Ejv„ F NnC^nl^ t \0) we obtain F Nn = for all N + 
0. Thus we have proved that 

|G> = G ^^J0> = |SVBS) . (C15) 

In summary, the state |SVBS) in Eq.(|4|) is the unique 
ground state of the generalized pseudo-potential Hamil- 
tonian ([1061) . 
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